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O ■ Abstract 

We simulate 3-flavour lattice QCD at finite temperature and isospin chemical 
Q-f potential In particular we study the nature of the finite temperature tran- 

^ ■ sition for quark masses close to the critical mass at which this transition for 

^ . zero chemical potentials changes from a first order transition to a crossover, 

^vq | We find that the Binder cumulants, used to determine the position of this tran- 

sition, have very strong dt dependence. This leads us to an estimate of the 
critical mass which is about 20% below previous estimates. In addition, when 
I/-) this dt dependence is taken into account, we find that the Binder cumulants 

Q\ • show very little dependence on fij. From this we conclude that we do not as 

■ yet see any evidence for the expected critical endpoint. We have argued previ- 

O '. ously that the position and nature of the finite temperature transition for small 

JO ■ fij should be the same as that for small quark-number chemical potential \i. 

1 Introduction 

p QCD at finite baryon-number density and temperature was present in the early universe and is produced 
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in relativistic heavy-ion colliders. QCD at finite baryon-number density and zero temperature describes 
nuclear matter, such as is found in neutron stars. 

Because QCD at a finite chemical potential /i for quark number has a complex determinant, stan- 
dard simulation methods, based on importance sampling, fail. At small /x, in the neighbourhood of the 
finite temperature transition, various approaches have been devised to circumvent these problems. These 
methods include reweighting [1], analytic continuation EH3JH and series expansions j5j|6]]. We adopt 
a simpler approach and simulate with at a finite isospin chemical potential m Q. This corresponds 
to simulating with finite quark-number chemical potential fi = /ir/2 using only the magnitude of the 
fermion determinant and ignoring its phase. For /i sufficiently small, this phase is well enough behaved 
on the sizes of lattice needed for the simulations, that one might expect that the position and nature of 
the transitions at finite /j and those at finite fij should be the same. 

For 3 flavours of quarks, at zero chemical potentials, there is a critical mass m c such that the finite 
temperature transition is first order for m < m c , second order in the 3-d Ising model universality class 
at m = m c and becomes a crossover for m > m c |8 ] (see also |9JI3)- Here it was expected that, as |/i| 
or \muj\ was increased, m c would also increase, providing a critical endpoint at some small chemical 
potential. 

We are performing simulations with m close to m c on 8 3 x 4 and 12 3 x 4 lattices. What we 
find is that the fourth order Binder cumulants, which are used to determine the nature of the transition, 
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are very sensitive to the updating increment dt used in the hybrid molecular-dynamics (R algorithm) 
simulations 1101 When this dependence is taken into account, m c ^ 0.027 compared with previous 
estimates of m c « 0.033. In addition we have determined that the Binder cumulants have very little /in- 
dependence. Thus, as yet, we have obtained no evidence for a critical endpoint. Both these observations 
are in agreement with the new results of de Forcrand and Philipsen flUl . 

In section 2 we describe our simulations and presents our results. Discussions and conclusions are 
presented in section 3. 



2 Simulations and Results 
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Fig. 1: Mass dependence of B^tyi/)) for various values of the updating increment dt. 



The staggered quark action for lattice QCD at finite \xj is 

Sf = Yj Ww^^i) + m]x + i\exT 2 x\ • (1) 

sites ~ 

which yields a real positive fermion determinant. We simulate this theory with 3 flavours of light dy- 
namical quarks on 8 3 x 4 and 12 3 x 4 lattices using hybrid molecular-dynamics (R algorithm) methods, 
in the viscinity of the finite temperature transition. Since we are interested in small /// (/i/ < m n ) we set 
A = 0. 

We run at several quark masses in the range 0.025 < m < 0.04, that is, close to the critical mass. 
The position of the transition is determined by the minimum of the fourth-order Binder cumulant [ 12 1 
for the chiral condensate, defined by 

B A = ^== K Z-" 1 , (2) 

where the overline indicates an average over the lattice. 5 noisy estimators of ipip are made at the end 
of each trajectory, which yields an unbiased estimator for B^ipip). We have checked that this method 
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Fig. 2: fii dependence of Bi^iptp) for various values of the updating increment dt at m = 0.035. 



for determining the position of the transition produces an estimate which is consistent with the position 
of the maximum of the corresponding susceptibility. This measurement is obtained by simulating at 
typically 4 values of j3 = 6/g 2 in the neighbourhood of the transition for a given mass, \ij (and dt), and 
continuing to the minimum using Ferrenberg-Swendsen reweighting [ 13 1. These simulations are usually 
160,000 length-1 trajectories long for each (m, /3, fii). 

Because B4 depends strongly on the updating increment dt, we have performed simulations for 
several different dts in the range 0.01 < dt < 0.0625 for each (m,^i). Figure ^ shows the Binder 
cumulants as functions of mass for the various dt values used, for our 12 3 x 4 runs. As we can see, the 
dt dependence is considerable. At the critical mass, B4 should take the Ising value (1.604(1)). If we 
estimate this from our dt = 0.05 simulations we would conclude that m c 0.0335. If on the other hand 
we use our dt = 0.02 runs we predict m c « 0.027. Thus we conclude that m c < 0.027. 

Figure[2]shows the [ii dependence of for m = 0.035 over the range < /xj < 0.3 (m^ ~ 0.4), 
also on a 12 3 x 4 lattice. What we see is that, for larger dt values, B4 appears to decrease with increasing 
[ii, which would lead us to the erroneous conclusion that there was a critical endpoint where B4 passes 
through the Ising value for some < 0.2. As dt — > 0, this falloff becomes less pronounced, and 
there is even the possibility that B4 increases with increasing /i/. Thus we find no evidence for a critical 
endpoint. 

In figure |3j we plot the values of B4 at m = 0.035 and [ii = against dt 2 . From this it is clear 
that, at the smallest dts, B^ is close to its dt = value. As we go to larger dt, B4 decreases making 
the transition appear more abrupt, i.e. less like a crossover and more like a first-order transition. Also 
included in this figure are the chiral susceptibilities, showing their strong dt dependence. These increase 
with increasing dt, again making the transition appear more abrupt. 

The transition 0, C and hence the transition temperature decrease slowly as fii is increased. Our 
best fit to our measurements of the fij dependence of f3 c at m = 0.035 from our dt = 0.02 runs gives 



(3 C « 5.15193 - 0.1758//?. 



(3) 
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Fig. 3: dt dependence of Bi^ip) and at m = 0.035 and fii = 0. The crosses are the Binder cumulants the circles are 
the susceptibilities. 



This is in reasonable agreement with the estimate at finite fj, of de Forcrand and Philipsen if we 

make the suggested replacement fij = 2fi. 



3 Discussion and Conclusions 

We are performing simulations of 3-flavour lattice QCD, using the hybrid molecular-dynamics (R) algo- 
rithm, at small isospin chemical potential, close to the finite temperature transition from hadronic/nuclear 
matter to a quark-gluon plasma. The quark mass is chosen to be close to the critical value for zero chem- 
ical potentials. The f3 and hence temperature of this transition decrease slowly with increasing in a 
manner consistent with the decrease with increasing fi, provided one identifies /i/ = 2/i. 

The nature of the transition is determined using Binder cumulants for the chiral condensate. We 
find that the these cumulants depend strongly on the updating increment dt, decreasing with increasing 
dt. This can be understood from the fact that the shift in effective (3 at finite dt is much larger below 
the transition than above it. Hence a small change in j3, which takes the system through the transition, 
produces a much larger change in the effective f3. This in turn induces a larger change in observables, 
which makes the transition appear more abrupt, which is reflected in a smaller value for B±. When this 
dependence is taken into account, the critical mass is found to be some 20% below previously published 
values. This is in agreement with the new results obtained by de Forcrand and Philipsen using RHMC 
(exact) simulations ATI . 

When we apply the finite dt corrections to our measurements of the Binder cumulants at finite /if, 
we find that they show little dependence on fij. So far, we find no evidence for the predicted critical 
endpoint. This too is in agreement with the work of de Forcrand and Philipsen, whose new predictions 
of the nature of the transition at finite \i obtained from continuations from imaginary also show a weak 
dependence of on /j \ 11 1. In fact, their predictions suggest that B^ might actually increase with /j as 
was suggested by some of our early results. 

We are continuing our simulations at a quark mass of m = 0.03 which is closer to the newly 
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determined critical mass, and at m = 0.025 which appears to be just below the critical value. We too are 
converting to the new RHMC algorithm [14], which is exact. That is, it produces results which do not 
have any finite dt errors. 

If no critical endpoint is found, it will indicate that the critical endpoint, if it exists, is unrelated to 
that found at zero chemical potentials, as the mass is varied. 

We intend to use simulations at finite [ii as a platform for reweighting (by the fermion phase 
factor) to finite \i. This has been found to be a better choice than reweighting from zero /i. Presumably 
this is because the position of the finite temperature phase transition is the same (or at least close) for 
the 2 theories. In this respect it has similarities with the method of reweighting from shifted (3 values, 
used by Fodor and Katz. Here we expect that new methods for simulating fermions will give us at 
least stochastic estimators for (the phase of) the fermion determinant, removing the burden of exact 
determinant calculations, which are too expensive. 
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